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ABSTRACT 

Narrow Line Seyfert 1 (NLS1) galaxies are remarkable for their extreme continuum and emission line 
properties which are not well understood. New results bearing on the spectroscopic characteristics of 
these objects are presented here, with the aim of establishing their typical ultraviolet (UV) and optical 
spectral behavior. We employ Hubble Space Telescope (HST) observations of 22 NLSls, which represent 
a substantial improvement over previous work in terms of data quality and sample size. High signal-to- 
noise (S/N) NLS1 composite spectra are constructed, allowing accurate measurements of the continuum 
shape and the strengths, ratios, and widths for lines, including weak features which are barely identifiable 
in other Active Galactic Nuclei (AGN) composites. We find that the NLS1 sources have redder UV-blue 
continua than those typically measured in other quasars and Seyferts. Objects with UV line absorption 
show redder spectra, suggesting that dust is important in modifying the continuum shapes. The data 
also permit a detailed investigation of the previously proposed link between NLSls and z > 4 quasars. 
Direct comparison of their composite spectra, as well as a Principal Component Analysis, suggest that 
high-z QSOs do not show a strong preference toward NLS1 behavior. 



Subject headings: galaxies: Seyfert — galaxies: abundances- 
lines 



galaxies: evolution — quasars: emission 



1. INTRODUCTION 

Narrow-Line Seyfert 1 galaxies (Osterbrock & Pogge 
1985; Goodrich 1989) are a subclass of AGNs that mani- 
fest a distinctive ensemble of properties. They are rather 
rare objects that exhibit relatively narrow broad lines, 
strong Fe II emission, and weak emission from the nar- 
row line region; they are more variable in X-rays than the 
Broad-Line Seyfert 1 objects ("normal" Seyfert Is, here- 
after Syl), and exhibit pronounced soft X-ray excesses. 
They seem to cluster at one extreme end of the Boroson 
& Green (1992) "Eigenvector 1" (EV1) relation, as a re- 
sult of their tendency toward weak [O HI] AA4959, 5007 
emission, and narrow and blue-asymmetric H/3 profiles. 
Understanding this EV1 is important for NLSls in partic- 
ular and for AGN in general since it may be closely linked 
and possibly driven by the central engine parameters, in 
particular L/ LEdd-, the Eddington ratio (Boroson & Green 
1992; Boroson 2002). However, whether this is the main 
and the only physical parameter that controls the NLS1 
classification, with its distinctive features, is still a matter 
of debate. 

To date, NLS1 studies are built on either individual ob- 
jects or on samples for which the spectra span only narrow 
wavelength ranges; also, with the exception of the X-ray 
observations, these samples are small in general. In par- 
ticular, detailed investigations of the NLS1 blue/UV emis- 
sion properties have been limited (e.g., Rodrigues-Pascual 
et al. 1997; Kuraszkiewicz & Wilkes 2000). Studies of 
their spectral energy distributions (SEDs) that cover wide 
bandpasses at a single epoch are almost completely miss- 
ing. This kind of data is particularly useful for testing and 
constraining models proposed for these sources. Further 
examination of larger samples of NLS1 emission spectra is 
clearly desirable. 



Understanding the nature of NLSls requires a detailed 
description of their average behavior. Therefore, a first 
goal of the present study is to obtain a comprehensive 
spectral characterization of the typical NLS1 galaxy. The 
mounting number of high quality HST spectra of NLS1 
sources allows for a better definition of their spectral prop- 
erties in general, and their UV emission in particular. In 
this study, we make use of HST archival observations of 22 
NLSls, a sample which is nearly twice as large as any of 
those used in previous studies of the NLS1 ultraviolet line 
and continuum emission. This database includes several 
objects whose observations cover a wide wavelength range 
(Lya to Ha region) permitting thus, for the first time, a 
simultaneous survey of the NLS1 UV and optical spectral 
features. We construct average and median NLS1 com- 
posite spectra and provide measurements of the resulting 
underlying continuum, and the line strengths, ratios, and 
widths. 

The second important question we attempt to address 
here is the degree to which NLSls and high redshift (z > 4) 
quasars share common properties. The possible connec- 
tion between these two classes of objects has been sug- 
gested by Mathur (2000), who proposed a scenario in 
which both NLSls and high z QSOs are in an early evolu- 
tionary phase, such that accretion proceeds at or near the 
Eddington limit. This analogy is mostly based on the indi- 
cations for high metallicities in both of these categories of 
sources, and on the presence of an enhanced low-velocity 
component in the UV spectra of high redshift quasars. 
However, no direct comparison of their emission properties 
has been attempted yet. This is primarily due to their wide 
separation in redshift, which makes it difficult to observe 
them in the same spectral ranges. The current availability 
of high-quality UV spectra for the low redshift objects al- 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space Telescope Science 
Institute. STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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lows us to test the validity of this picture. Wc examine and 
compare the emission properties of the NLSls and z > 4 
quasars, in order to determine the extent of the similarities 
exhibited by these two classes of objects. This comparison 
provides additional insights into AGN behavior as a func- 
tion of redshift, luminosity, metallicity, and other physical 
parameters. 

2. THE NLSl SAMPLE & DATA PROCESSING 

Non-proprietary, archival HST Faint Object Spectro- 
graph (FOS), Goddard High Resolution Spectrograph 
(GHRS), and Space Telescope Imaging Spectrograph 
(STIS) spectra of all objects known to us that were pre- 
viously identified in the literature as NLSl galaxies (i.e., 
FWHM(H/3) < 2000 km s" 1 , [O IIl]/H/3 < 3), were re- 
trieved from the Space Telescope Science Institute in the 
form of calibrated data. A number of these objects are 
borderline sources in terms of this specific classification; 
however, they resemble the NLSl characteristics by ev- 
ery other definition, and they have been extensively used 
in other NLSl studies. Table 1 summarizes the instru- 
mental setup and the resulting total wavelength coverage 
corresponding to each observation. References relevant to 
the classification of these objects are also listed. Each line 
in the table refers to observations obtained under a sin- 
gle observing program; in all but one case (I Zw 1), these 
observation sets were taken at a single epoch. The data 
obtained under separate observing programs are listed as 
different lines. 

In general, observations of each object were acquired 
with multiple gratings. The full wavelength coverage was 
obtained by co-adding the individual spectra, after re- 
sampling to a common dispersion (the lowest number of 
A/pixel, to avoid loss of information), and flux-averaging 
in the overlapping regions. As indicated in Table 1, there 
are several galaxies for which multiple observations were 
obtained at different epochs and in some cases with dif- 
ferent instruments. Because these sources are character- 
ized by significant variability, the continuum level often 
needed minor rescalings (< 10%) before averaging. The 
use of multiple instrumental configurations translated also 
into different spectral resolutions R. To obtain the final 
observed spectrum, it was necessary to convert the individ- 
ual observations first to a common R, by gaussian smooth- 
ing of the spectra with high resolution, and second, to a 
common dispersion. For all objects, the averaging in the 
common wavelength ranges was performed by weighting 
the spectra by the reciprocal of their noise variance. 

Many of these spectra display significant resonance line 
absorption. The main focus of this analysis is on the emis- 
sion lines, and therefore, accurate measurements of these 
features required correction for the intervening absorption. 
When this appears as narrow lines superimposed on cither 
the emission lines or the continuum, we removed the ab- 
sorption feature through a simple interpolation. When 
the absorption was severe, especially near or within some 
of the strong emission features (Lya, N V, and C IV), a 
conservative reconstruction of the line was attempted by 
low-order polynomial fitting. For these particular cases, 
as indicated in Table 2, the line flux and equivalent width 
(EW) measurements are considered to be only lower limits. 

The objects in this sample are subject to small amounts 



of Galactic foreground extinction, spanning Ay = 0.03 — 
0.31, with a typical value of ~ 0.05 (Schlegel, Finkbeiner 
& Davis 1998). We corrected the spectra for the resulting 
reddening using the empirical selective extinction function 
of Cardelli, Clayton & Mathis (1989). No correction for 
intrinsic reddening was attempted. 

In order to combine the individual source observations 
into a single composite spectrum, a proper alignment in 
wavelength of the emission lines is necessary. This re- 
quires, in turn, consistent measurements of the redshift 
values z used in the Doppler-correction process. There- 
fore, we remeasured the redshifts using the C IV A1549 
emission-line, which is measured in 21 out of 22 objects 
(the C IV feature is a doublet AA1548.2, 1550.7, with the 
simple average of A1549.5, corresponding to the optically 
thick case, consistent with photoionization models). The 
observed wavelength for each case was established by fit- 
ting a Gaussian to the top 20% of the profile. For situa- 
tions where the measured line was significantly affected by 
absorption or low S/N, we verified the result with fits to 
the top 50% of the profile. The resulting z determinations 
are in good agreement with previously published values, 
with an estimated uncertainty of < 0.001, and lead to a 
good superposition of the principal emission lines in the 
rest-frame. The redshift measurements obtained for each 
object from C IV line fitting are listed in Table 2. 

Measurements of the emission features in the individ- 
ual spectra were carried out with line profile fitting. This 
was in general unsuccessful when single Gaussians were 
used, as they tend to lose flux from the wings if they are 
prominent, and from the peak if the cores dominate. The 
overall shapes of the lines measured in these NLSl spectra 
were generally well represented by single Lorentzian fits, 
which is consistent with previous findings (Moran, Halpern 
& Helfand 1996; Leighly 1999; Veron-Cetty, Veron & 
Goncalves 2001). The C IV equivalent widths are also 
listed in Table 2, along with the rest- frame continuum lu- 
minosity measured at 1450A [£,,(1450), expressed in ergs 
s -1 Hz -1 ]. The EW measurements may have errors result- 
ing from the choice of continuum placement; experiments 
with alternative continuum fits suggest a systematic un- 
certainty of ~ 15%. 

Another parameter recorded in Table 2 is the spectral 
index a of the power-law fit (defined by F v oc v a ) that 
best approximates the continuum shape of each individual 
object spectrum in the rest-frame. In determining the con- 
tinuum solution, we tried to make use of the wavelength 
ranges which contain pure continuum emission. A reliable 
fit is obtained when a wide separation of the continuum 
windows is available (e.g., ~ 1100 - 4000 A), such that the 
Fe emission, which is strong and prevalent in the UV and 
the blue part of the spectra of NLSls, is avoided. In the 
present sample, there is only a small number of objects 
for which the available spectra cover completely and/or 
extend redward of the broad Fe II and Fe III emission-line 
complexes. For the sources spanning only the UV range, 
the measured spectral index should be considered a lower 
limit, as the apparent continuum can be strongly contam- 
inated and reddened by the Fe emission. The full range of 
slope values and their distribution among the sample are 
displayed in Figure 1, along with the calculated median 
and average values. In identifying the Fe features, and con- 
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sequently the emissionless continuum windows, we used 
information based on high resolution spectra of I Zwl, as 
provided by Vestergaard & Wilkes (2001) in the UV band- 
pass, and by Oke & Lauer (1979) in the regions redward 
of Mg II A2797.9. The windows chosen in determining the 
power law continuum shape are placed around 1140, 1285, 
1320, 1350, 1450, 3810, 3910, 4040, 4150, 5470, 5770-5800, 
and 6210A, with the interval lengths ranging from 10 to 
20A. These wavelength ranges are also used in fitting the 
composite spectra (see Section 3.1), and they are illus- 
trated in Figure 2. 

3. THE NLSl COMPOSITE SPECTRA 

In this section we present NLSl composite spectra, and 
resulting measurements of continuum and emission-line 
characteristics. The use of composite spectra is comple- 
mentary to analysis of individual spectra. Composites of- 
fer higher S/N ratios, allowing measurement of weak fea- 
tures, while also directly providing a description of typical 
NLSl spectral properties. We examine the results in re- 
lation to similar observations of broader AGN samples, in 
order to explore the extent to which NLSls represent a 
distinct subclass. These comparisons potentially provide 
useful tests for discriminating between model scenarios for 
the NLSl phenomenon. 

3.1. Overall continuum 

When building spectral composites, one of the most 
challenging tasks is the generation of an underlying con- 
tinuum that reflects the typical appearance of the sample 
as a whole. Francis et al. (1991) noted that "there is no 
'correct' way of co-adding spectra that exhibit differences 
on many different scales." Our sample consists only of 
NLSl sources, which, as a subcategory of the AGN clas- 
sification, share a set of common attributes. The spectra 
nonetheless show a substantial diversity in their proper- 
ties. As Figure 1 and Table 2 reveal, the continuum shape, 
as described by the power-law index a, exhibits a consider- 
able variety among individual NLSl spectra. Because only 
very few objects in the sample span the whole wavelength 
range, a simple median or average composite displays ar- 
tificial discontinuities at locations where the number of 
contributing spectra changes significantly. Some care is 
thus required in order to use these composites to describe 
typical continuum properties for NLSl galaxies. 

Figure 2 presents the NLSl average and median compos- 
ites, constructed using data from individual objects nor- 
malized to the mean flux in the wavelength range 1430 A- 
1470 A. Each spectrum was given equal weight, thus avoid- 
ing biasing the resulting composites toward the brightest 
objects (i.e., those with the highest S/N). The median dis- 
plays a smoother continuum than does the average, and is 
well described by a single power law (a = —0.798 ± 0.007) 
from just redward of Lya to H7. Owing to the relatively 
high S/N of the spectrum and the wide separation of the 
fitted regions, the statistical uncertainty in the spectral 
index, as given by a Chi-square minimization method, is 



quite small. However, the value itself is rather sensitive 
to the precise wavelength sections employed in the fitting. 
Redward of H7, the continuum flux density rises above the 
level predicted by the UV power-law, and is best approx- 
imated by a separate power-law, with a = —2.38 ± 0.01. 
Since these objects are low-luminosity AGNs, contamina- 
tion by the host galaxy starlight may contribute to this 
change in the spectral shape. Both fits are shown in the 
middle panel of Figure 2; the continuum windows used in 
the power-law fits are indicated as horizontal lines below 
the composite. Prominent emission features in the opti- 
cal range are labeled. The number of objects contributing 
to the composites at each wavelength is presented in the 
bottom panel of Figure 2. 

The UV-bluc continuum spectral index for this HST 
NLSl composite falls among the steepest values found in 
other AGN composite measurements (—1 < a < —0.4; 
see for example, Table 5 in Vanden Berk et al. 2001) 2 . 
Such red continua were measured only in the Zheng et 
al. (1997) and Telfer et al. (2002) HST composite spec- 
tra (a = —0.99, and a = —0.71 respectively), and in the 
Schneider et al. (2001) sample of Sloan Digital Sky Survey 
very high redshift quasars (average a = —0.93). There are 
several possible reasons why these particular samples show 
such steep continua, and we explore them in the following 
paragraphs, in an effort to understand the origins of the 
extreme continuum properties of the NLSl composites. 

One effect that may contribute to the soft continuum 
slopes of the high z QSOs and the HST sources derives 
from the restricted wavelength baseline typically used in 
fitting the continua. The spectra of objects comprising 
these steep continuum samples make available only short 
regions redward of Lya, in which the iron contamination 
may produce an artificial rise in the continuum profile 
(Telfer et al. 2002) . However, in the NLSl composites, 
the wide baseline used in fitting the continuum permits a 
successful accounting for the Fe II and Fe III complexes, 
therefore making the Fe contamination an unlikely cause 
for their overall red continua 3 . 

The steepness of the HST composites can also poten- 
tially be attributed to an evolution of QSOs to softer spec- 
tra at lower redshifts and/or an observational bias toward 
detecting sources with harder continua at higher redshift 
(Francis 1993). The space-based observations include a 
significant number of low redshift sources, and this is the 
case for the NLSl galaxies as well. The same evolution- 
ary effect and the potential detection biases can also be 
responsible, at least partially, for the redness of the NLSl 
composites, especially when compared with the ground- 
based composite spectra. 

The present sample shows little evidence of a correla- 
tion between continuum slope and z, but a considerably 
stronger trend relating a and luminosity (Fig. 3). Regard- 
less of its origin, the luminosity dependence of spectral 
slope is probably largely responsible for the relatively steep 
a for the NLSl composite, since the typical luminosity of 
our NLSls is lower than that of sources employed in most 
other quasar composites. 



2 Rodrigucs-Pascual ct al. (2000) have also reported evidence that NLSls are redder than other Seyfert Is at optical wavelengths. 

3 The Fe contamination is the primary factor that accounts for the difference between the spectral index that best fits the NLSl median com- 
posite and the median (and average) value of the power- law index distribution of the individual continuum fits (see Figure 1). The majority 
of the objects in the NLSl sample cover only the UV wavelength region, where the power-law solution is unable to properly estimate the true 
underlying level below the Fe II and Fe III emissions, and as a consequence, steeper indices are measured. 
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3.2. Reddening and Absorption 

The luminosity dependence of spectral slope could be 
intrinsic to the accretion source, but independent evi- 
dence suggests that it is mostly attributable to luminosity- 
dependent reddening. Internal dust, if present, is expected 
to be accompanied by gas producing observable absorption 
lines. In our sample, the signature of strong absorption 
near the systemic velocity of the host galaxy is present 
in almost half the objects (see Table 2). Median val- 
ues of the spectral indices corresponding to subsamples of 
sources with and without strong resonance line absorption 
are —1.34 and —0.73 respectively, consistent with steeper 
(redder) continua in the absorbed NLSls. Median loga- 
rithmic luminosity values for the absorbed and unabsorbed 
subsamples are 29.01 and 29.60 respectively (mean values 
are 28.4 ± 0.4 and 29.3 ± 0.2 respectively), thus directly 
linking the presence of absorption with luminosity. 

The location and physical state of the absorbers in 
NLSls may be important for understanding these objects. 
Intrinsic absorption related to the central accretion source 
is now known to be present in a significant fraction of 
Seyfcrt 1 nuclei (Crenshaw et al. 1999). NLSls exhibit 
similarities in emission properties to low-ionization broad 
absorption-line QSOs (BALQSOs), which Boroson (2002) 
has interpreted as evidence for high luminosities relative 
to the Eddington value for both classes of object; in the 
BALQSOs, the absorption is clearly also closely related 
to the accretion process. In the NLSls, the intervening 
gas may be a warm (highly ionized) absorber that is po- 
tentially dusty (e.g., Komossa & Greiner 1999); alterna- 
tively, the absorbing matter may be related to the "luke- 
warm absorber" identified in case studies by Kraemer et 
al. (2000) and Crenshaw et al. (2002), which may reside 
on kpc scales. The luminosity dependence of reddening 
and absorption for the NLSls is consistent with a sce- 
nario in which the absorbing matter covers a larger solid 
angle as seen by the central source in lower luminosity 
objects, a natural expectation for interstellar matter in a 
disk- like geometry of luminosity-independent scale-height, 
but with a sublimation radius for dust that scales with lu- 
minosity, similar to the 'receding torus' model (Lawrence 
1991). If the absorbing medium has a flattened distri- 
bution aligned with the host galaxy disk, we might then 
expect a correlation between the host inclination and the 
spectral (UV) color, as reported for several Seyfert Is by 
Crenshaw & Kraemer (2001). Inclination values are avail- 
able in the literature for 16 of our sample members; the 
data do not show a statistically significant trend with a, 
however. In summary, the extent to which the UV/optical 
absorbers in NLS1 nuclei are associated with the central 
accretion structure versus the normal host galaxy inter- 
stellar medium remains ambiguous. 

The nature of the dusty absorber in these galaxies can 
be further tested by analyzing their soft X-ray (0.1 - 2.4 
kcV) characteristics. A "lukewarm" (low ionization) con- 
stitution of the gas would suggest strong absorption, that 
translates into flatter observed soft X-ray spectra. If a 
"warm" (high ionization) absorber is present, steeper X- 
ray continua should be measured (Grupe et al. 1998). 
The R OS A T spectral indices (L = 1 — a) are available for 
the majority of the NLSls employed in this study, and are 
recorded in Table 2. Splitting the sample into absorbed 



and unabsorbed objects, as before, and calculating the me- 
dian value of the X-ray spectral slope for each subgroup, 
should provide us a rough criterion for distinguishing be- 
tween the two types of absorbers. The resulting values are 

Tabsorbed spectra — 3.4 and F una bsorbed spectra 3.1, nomi- 
nally implicating the presence of highly ionized absorbers 
(mean values are 4.2±0.8 and 3.3±0.2 respectively); how- 
ever, the difference in the median values is small, and it 
appears likely that the absorbing medium in these objects 
is described by a range of properties. 

3.3. The UV spectrum 

Figure 4 shows the average spectrum in the UV range, 
where most of the objects (> 18) in the sample are con- 
tributing. The standard deviation (RMS) and the stan- 
dard deviation of the mean spectra show the degree of 
variation between the members of the sample and the un- 
certainty of the average composite as a function of wave- 
length. Near the wavelength corresponding to the common 
continuum normalization (A ~ 1450A), the spectral vari- 
ation within the NLS1 sample is dominated by the mod- 
ulations in the cores of the strongest emission lines (Lycn, 
C IV, He II, and C III]). In this respect, NLSls behave 
like the other more general samples of AGNs (e.g., Fran- 
cis et al. 1992; Brotherton et al. 2001). Away from the 
normalization wavelength interval, the bulk of the spectral 
variance is accounted for by differences in the individual 
continuum shapes and noise. 

3.4. The optical composite 

For wavelengths A > 3000A, the fiSTNLSl spectrum is 
based on only a small number of objects. Because fewer 
than 5 NLSls make up the composites at these wave- 
lengths, the individual source contributions are more pro- 
nounced, making it difficult to identify and measure the 
spectral features (continuum and emission lines) located 
at the transition points in the number of contributing 
objects (e.g., at ~ 4360A, ~ 5500 A). A smoother spec- 
trum is constructed using only the three spectra that span 
the whole wavelength range (Ark564, Mrk493, WPVS007), 
and is presented in Figure 5. We compare our optical me- 
dian composite with the Sulentic et al. (2002) median 
NLS1 spectrum, constructed from a much bigger sample 
(24 sources, ground-based spectra). The two spectra are 
very similar in both the continuum shape and the pro- 
file and strength of emission features. Small differences, 
like the stronger peaks in the emission lines in the HST 
composite, are mostly due to the different spectral resolu- 
tion that characterizes the two samples ( 4 - 7A FWHM 
for the Sulentic spectra, ~ 2A for the FOS-iJST obser- 
vations). The good agreement between the two medians 
suggests that the HST composites are representative of 
typical NLS1 optical spectra, despite the fact that they 
are constructed from such a small object sample. 

3.5. Emission Lines 

The richness of the line emission in the AGN spectra 
is easily distinguishable in the NLSls due to their narrow 
widths. Hence, careful identification and intensive analysis 
of the NLS1 emission lines have been performed, though 
to date, only on individual sources, e.g. I Zwl (Laor et 
al. 1997; Vestergaard & Wilkes 2001, and references 
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therein) . When composites are built using relatively large 
samples, higher S/N is achieved, and therefore, higher ac- 
curacy is expected in identifying and measuring the emis- 
sion lines. We have been able to detect and parametrize 
in the HST NLS1 composites a large number of emission 
features. Many are barely present or are heavily blended, 
and therefore difficult to quantify in other AGN compos- 
ites (Francis et al. 1991; Vanden Berk et al. 2001; Zheng 
et al. 1997; Telfer et al. 2002). 

The emission lines are measured from the average com- 
posite, as it presents a higher S/N ratio than the median 
spectrum. An exception is made for the wavelength in- 
terval ~ 2900 — 3170 A, where the emission features in 
the average spectrum are distorted by the jumps present 
in the continuum shape (see Section 3.1), and where the 
median is used. Line fluxes, strengths (EWs), widths 
(FWHM) and line shifts relative to the laboratory cen- 
tral wavelengths (Av) are listed for all detected features, 
along with their 1 a error bars, in Table 3. The lines have 
been identified by matching wavelength positions and re- 
spective relative strengths of features found in the Francis 
et al. (1991), Zheng et al. (1997), Vanden Berk et al. 
(2001) and Telfer et al. (2002) composites, and in the 
detailed analyses of I Zwl (Laor et al. 1997; Vestergaard 
& Wilkes 2001). Tentative identification of newly discov- 
ered features is based mainly on data available from the 
Atomic Line List 4 and Verner et al. (1996). 

All measurements, including the peak positions of each 
emission line, A mea „_ rest , and the respective errors are gen- 
erated using the task specfit (Kriss 1994) as implemented 
in the IRAF 5 software package. The method employs line 
and continuum spectral fitting via an interactive x 2 mini- 
mization. The fit is performed in 14 separate spectral in- 
tervals; their lengths varied, in order to keep a reasonable 
number of parameters, from ~ 150A in the UV-blue re- 
gion, rich in prominent emission features, to ~ 500 - 900A 
in the optical range. The equivalent widths are measured 
relative to the resulting local continua, which differ from 
the shape estimated in Section 3.1, mainly due to the fact 
that the average and not the median composite is used in 
the fitting process. Comparisons indicate that line fluxes 
measured in the average and median composites are con- 
sistent to within < 10% in most cases. 

For most features, the line emission is modeled using 
single Lorentzian profiles. However, the most prominent 
lines required two components, a narrow and a broad 
Lorentzian, the broad one being in general blueshifted rel- 
ative to the narrow component. The broad and prominent 
Fe emission surrounding and redward of Mg IIA2798 is 
not measured; thus, line measurements for 2000 A < A < 
4250A may include some Fe contamination, and therefore, 
a high level of uncertainty. For the optical wavelength 
range, redward of H7, we use in the fitting process the 
empirical Fe II template obtained by Boroson & Green 
(1992) from the I Zwl spectrum. No broadening of the 
iron template was necessary in order to match the line 

4 The Atomic Line List is hosted by the Department of Physics and Astronomy at the University of Kentucky; 
http: / /www. pa. uky.edu/~peter/atomic. 

5 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy Inc. (AURA), under cooperative agreement with the National Science Foundation. 

6 Single Lorentzian fits for these features result in line shifts that differ by ^100 km s _1 from the ones given by the narrow component only; 
this is due to the fact that the cores of these relatively narrow lines account for the bulk of the flux. 

7 The velocity shifts may be unusually large in I Zw 1; see Laor (2000). 



width of the NLS1 composite spectra. In fitting forbidden 
line doublets, the line pairs are assigned to have common 
velocity widths and offsets, and their flux ratios are con- 
strained to values determined by branching ratios, when 
appropriate. 

Because the composites were constructed using redshifts 
based on the position of a single emission-line, C IV (see 
Section 2), a check for systematic velocity offsets for other 
lines can be performed. Such line shifts, which are def- 
initely present for many of the lines listed in Table 3, 
have been detected previously in other AGN samples (e.g., 
Gaskell 1982; Espey et al. 1989; Tytler & Fan 1992; 
Mcintosh et al. 1999; Vanden Berk et al. 2001). Note 
that most of the lines are redshifted relative to C IV, the 
reference line, and that the largest values of the recorded 
velocity shifts correspond to the lowest ionization states 
(e.g., Si II, C II, O I). Figure 6 shows the velocity shifts 
versus the ionization potentials, compiled for all emission 
lines that are stronger than C IIA1335, i.e., EW > 4.0 A, 
and that have well defined non-blended peaks. We treat 
the permitted and the forbidden lines separately since they 
may have very different origins. The values for the ioniza- 
tion potentials are chosen such that, for the recombination 
lines of H and He, they express the energy necessary to 
ionize the respective state for later recombination, while 
for the collisionally excited lines from heavy elements, they 
represent the energy that is needed to create the ionization 
state. For the most prominent features, which are fitted 
with two components, we consider only the velocity shift 
measurements given by their narrow components (they are 
identified by different symbols in Figure 6, upper panel); 
this treatment is justified by the fact that the reference 
wavelength is chosen based on the location of the peak of 
C IV, which is measured using the narrow component fit 6 . 
This approach is consistent with measurements of velocity 
offsets performed in previous similar studies of other AGN 
samples (e.g. Vanden Berk et al. 2001). 

The results presented in Figure 6 provide evidence for 
an anticorrelation between the velocity shifts and the de- 
gree of ionization of the emitting species. The trend is 
definitely present in both categories of emission features. 
For the forbidden lines, significant velocity shifts, usually 
blueshifts relative to the systemic velocity, have been de- 
tected in the past for the high ionization transitions (e.g., 
Penston et al. 1984; Appenzeller & Wagner 1991); how- 
ever, the evidence for a correlation similar to that pre- 
sented in this study has only recently been revealed by 
measurements from the high S/N Sloan Digital Sky Survey 
quasar composite spectrum (Vanden Berk et al. 2001). In 
both types of lines measured in the NLS1 sources, the am- 
plitude of this correlation appears similar to that present 
in other, more heterogeneous Syl/quasar samples 7 . Al- 
though the ionization potential seems to govern the magni- 
tudes of the velocity offsets, there is a significant amount of 
scatter in this empirical relationship, in both Syl/quasars 
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and NLS1 objects, suggesting the influence of other pa- 
rameters as well. The scatter in the NLS1 trend may be 
amplified by wavelength zero-point calibration uncertain- 
ties, particularly for the FOS spectra. 

The origin of the velocity offsets between the AGN emis- 
sion lines is not well understood, but a promising means 
of interpreting this behavior is in terms of the disk-wind 
model for the broad- line region (BLR). In this picture, the 
high ionization lines are produced in outflowing winds, ac- 
celerated by radiative line driving, that arise from the ac- 
cretion disk, which is the base for the low-ionization emis- 
sion (e.g., Murray & Chiang 1998; Proga, Stone & Kali- 
man 2000; Lciglily 2001). Additionally, the disk and/or 
the radiative outflow itself, believed to proceed from near 
the plane of the disk, is assumed to be optically thick, so 
the emission from the receding wind is obscured. 

The fact that NLSls and other broad-line AGNs show 
similar line shifts provides additional constraints for the 
physical models that best characterize the NLSls. The 
preferred explanations for the primary drivers of the 
NLSl's extreme measured properties postulate either 
higher ratios of their luminosity to the Eddington lumi- 
nosity (e.g., Pounds et al. 1995), or a more pole-on 
view of an assumed disk-shaped BLR (e.g., Brandt & Gal- 
lagher 2000). The velocity offsets between lines reported 
here, and their similarity to those found in less restricted 
AGN samples, may present a challenge to either NLS1 sce- 
nario, if the velocity differences stem from a disk-wind phe- 
nomenon or something similar. For the orientation model, 
it would be surprising if the typical line-of-sight differences 
in line velocity remain essentially unchanged when the ori- 
entation is such that the line widths are diminished. For 
the high Eddington ratio picture, the offsets ultimately 
trace the disk outflow velocity, which is expected to scale 
with the local orbital speed, hence large velocity differ- 
ences combined with small widths are again unexpected. 
The solution in either case may ultimately involve the fact 
that the narrow line widths characterize only a subset of 
emission features in NLSls, notably including the H/3 fea- 
ture, while other lines (particularly in the UV) still re- 
veal high-velocity gas components. Indeed, as discussed 
by Wills et al. (2000), and elaborated by Shang et al. 
(2003), the UV line widths do not correlate in a simple 
way with H/3 width or with other NLS1 defining proper- 
ties. 

4. ARE NLSlS THE ANALOGUES OF HIGH Z QSOS? 

Based on the apparent similarity of some of the emission 
properties of NLSls and z > 4 QSOs, Mathur (2000) sug- 
gested a connection between these two classes of AGNs, in 
the sense that the NLSls are the low redshift, low lumi- 
nosity analogues of the high z quasars. In the proposed 
picture, both categories of sources are in an early evolu- 
tionary phase, in which accretion proceeds at or near the 
Eddington limit. This scenario has several appealing as- 
pects for explaining NLS1 phenomena, but additional tests 
are desirable to verify this idea, and especially to gauge its 
applicability to the high z sources. 

One of the key arguments employed in support of a 
NLS1 - high z QSO association was the initial report of 
narrow UV lines in z > 4 quasars by Shields & Hamann 
(1997). While the presence of an enhanced low- velocity 



component in the high redshift quasars was later confirmed 
by Constantin et al. (2002), the UV line profiles do not 
necessarily provide an adequate basis for linking the high 
z QSOs to a NLS1 classification, as this is based on optical 
lines. In fact, as noted above (Section 3.5), the UV and 
optical emission-line properties are substantially indepen- 
dent (Shang et al. 2003). 

Another possible similarity between high z quasars and 
NLS1 galaxies is a high metallicity (Z) in their emitting 
gas. In quasars at z > 3, super-solar heavy element en- 
richments have been derived in specific abundance studies 
based on both emission (Dietrich et al. 1999; Dietrich & 
Wilhclm-Erkens 2000; Dietrich et al. 2003) and absorp- 
tion properties (e.g., Hamann 1997). The high gas-phase 
metallicity is believed to derive from a recent episode of 
vigorous star formation. For the NLSls, several lines of 
argument have suggested the presence of enhanced metal- 
licity. One of the potential indicators of high abundances 
in these objects is their characteristically strong Fe II emis- 
sion (Collin & Joly 2000). However, the excitation of this 
emission is complicated and its strength does not necessar- 
ily map in a simple way to the Fe abundance (e.g., Verner 
et al. 1999). Moreover, if the Fe is largely produced in 
Type-la supernovae, as would be expected, the evolution 
of the progenitor population and the resulting enrichment 
should occur on a timescale > lGyr; consequently, a high 
iron enrichment is not a strong indicator of youth, and 
therefore, would not support the hypothesis of NLSls be- 
ing galaxies in the making. 

A better tracer of abundances in NLS1 galaxies is pro- 
vided by nitrogen. This clement is of particular interest 
since its production is believed to be dominated by sec- 
ondary enrichment, which translates into N/H oc Z 2 . Na- 
gao et al. (2002) have recently discussed the forbidden line 
spectra of NLSls, including [N II] A6583, and argued that 
the observed line ratios suggest higher metallicitics than in 
the Syls; as the authors also point out, however, this state- 
ment is very model-dependent. More robust diagnostics 
are obtained from studies of the N V A1240 feature, and 
in particular, measurements of its strength relative to C IV 
A1549 and He II A1640 (Hamann et al. 2002, and refer- 
ences therein) . Evidence for systematic trends of stronger 
N V and weaker C IV with increasing EV1 (which seems 
to trace NLS1 behavior, see Section 1) has been reported 
by Wills et al. (1999). Based on the same flux ratios, 
Shemmer & Netzer (2002) showed that the NLS1 sources 
deviate significantly from the well known relationship be- 
tween Z and luminosity (L) in AGNs (Hamann & Ferland 
1993), by exhibiting higher Z at a given L. Interestingly, 
several of the nine extreme NLSls used in their study have 
line ratios as high as those measured in the most luminous 
high z QSOs. Measurements for larger samples of NLSls 
are clearly desirable to verify these findings. 

In this section, we investigate the NLS1 - z > 4 QSO 
analogy by directly comparing their emission-line proper- 
ties. As most of the line emission studies of the NLS1 
objects are in the optical regions, and the corresponding 
(rest- frame) data for z> 4 are almost nonexistent, the UV 
spectral observations are the only accessible tool to ex- 
ploit in such a comparison. We have available both NLSls 
and high z quasar UV data, for samples that permit an 
adequate statistical analysis of their emission characteris- 
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tics. Our study employs both comparisons of their com- 
posite spectra, and an eigenvector analysis that determines 
the degree to which the spectral variances throughout the 
samples share common properties. The NLS1 objects and 
composites used here are described in Sections 2 and 3.1. 
For the z > 4 QSO sample, we made use of 44 spectra of 
non-BAL quasars, presented by Constantin et al. (2002) 
and Schneider, Schmidt & Gunn (1991), which span the 
1100 - 1700A rest-frame range. 

4.1. Direct comparison of composite spectra 

Figure 7 shows the NLS1 average composite (Sec- 
tion 3.1) ovcrplotted on the z > 4 QSO average compos- 
ite (Constantin et al. 2002). In a first approximation, 
the spectra agree well: similar continuum shape and the 
same strong emission lines, with comparable profiles. The 
pronounced Lya forest in the z > 4 composites does not 
reflect an intrinsic difference in the nature of the emission 
sources, but the expected increase in the opacity of the 
intergalactic medium at larger redshifts. 

The most evident distinction between the NLS1 and 
high z quasar composite spectra is in the strength of 
the principal emission features: Lya, Si IV+O IV], C IV 
and He II. This difference in the line strengths is a clear 
manifestation of the Baldwin Effect (the inverse equiva- 
lent width - luminosity relationship, Baldwin 1977). The 
two samples differ significantly in their average luminosity, 
with the lower luminosity objects, the NLSls, exhibiting 
stronger lines. Figure 8 shows the individual measure- 
ments of the rest-frame continuum luminosity, L„(1450A), 
and the strength (EW) of the C IV line, plotted for both 
the NLSls and high z quasars. The Baldwin Effect corre- 
lation in C IV is clearly present over the combined sample; 
the segregation in luminosity of the two samples is also 
evident 8 . Besides the basic Baldwin trend, the compar- 
ison in Figure 7 shows also the known trend with ion- 
ization: steeper Baldwin Effect for more highly ionized 
species (e.g., Dietrich et al. 2002). The difference in the 
strength of the emission lines in the NLS1 and high z QSO 
composites is likewise most obvious in Lya and the high 
ionization features (Si IV+O IV], C IV and He II), and 
almost absent in the low ionization lines (O I, C II). As 
found in previous studies of the Baldwin relationships, the 
N VA1240 emission proves to be the exception for which 
the line strength remains nearly independent of luminosity. 

Another important issue that can be examined via com- 
parison of the composite spectra is the degree of similarity 
in the chemical enrichment for the two source types. As 
noted before, the broad line region enrichments in both 
z > 4 QSOs and NLSls have been evaluated based on the 
N V/C IV and N V/He II line ratios, as best indicators of 
the overall mctallicity. Figure 7 shows that the N V emis- 
sion feature is equally strong in the NLS1 and high z QSO 
average composites. Deblended measurements of this line 
obtained from detailed fits of the Lya+N V+Si II*+Si II 
emission- line complex confirm this result. C IV and He II 
are however much stronger in the NLS1 spectra than in 
those of the z > 4 QSOs 9 . These trends result in lower 
N V/C IV and N V/He II line ratios, and presumably lower 

8 The flattening of the Baldwin correlation at low luminosities has bc< 
discussion and references. 

9 The same results are obtained when median composites are used for 



average metallicities, in the NLS1 sources than in the high 
z quasars, and thus do not support a strong connection be- 
tween the two classes of objects. The results of this com- 
parison differ from those of Shcmmer & Netzer (2002) who 
argued that the line strengths suggest comparable metal- 
licity for the NLSls and the high luminosity QSOs. How- 
ever, their conclusion was based on a very small number 
of NLSls; also, the high abundances in these sources were 
supported by a relatively large N V/C IV ratio, but not 
clearly confirmed by N V/He II. The discrepancy between 
their findings and ours may be related to the details of the 
samples used in each case. With the larger sample size em- 
ployed here, we tentatively conclude that the abundances 
in NLSls and high z QSOs are in fact characteristically 
different; further study with larger sets of object spectra 
would clearly be desirable. 

4.2. Principal Component Analysis 

Additional information on the NLS1 - z > 4 QSO re- 
lationship can be obtained by the means of a Principal 
Component Analysis (PC A), a mathematical decomposi- 
tion of a set of properties describing the sample into a 
smaller number of eigenvectors that can account for the 
bulk of the total variance present in the data. Due to the 
fact that quantitative measurements of the spectral prop- 
erties are in general subjective to the chosen parametriza- 
tion, we chose to apply the PCA method directly to the 
observed spectra (e.g., Francis et al. 1992). We proceed 
with this analysis for the two samples separately, and for 
the combined set as well. This approach should allow for 
further identification of potential common behaviors and 
their distribution among the two object classes. 

For the NLS1 data, the spectral range with the most ex- 
tensive coverage among the sample objects (21 out of 22) 
spans the region around the C IV line, and this bandpass 
is thus optimal for detailed statistical investigation. In 
this wavelength range, the z ^ 4 QSO sample is well rep- 
resented by all 44 object spectra. Therefore, we conduct 
the PC analysis on this interval, where the highest S/N 
is expected. Figure 9 shows the mean spectra, the RMS 
spectra, and the first five principal components (PCs), or- 
dered by the fraction of the sample variance for which each 
accounts, for the NLSls, high z QSOs, and the combined 
sample. The range of properties exhibited by these ob- 
jects is illustrated by the PCs and by their contribution 
to the total spectral variation. Principal Component 1 
(PCI) is dominated by line-core modulations (compared 
with the emission in the average composites), and there- 
fore, it can be considered primarily as a measure of the 
strength of the lines. This is a similar result to what Fran- 
cis et al. (1992) obtained with PC analysis for a larger and 
more heterogeneous sample (232 objects), from the Large 
Bright Quasar Survey (LBQS). Aside from a scale factor, 
this first principal component is very similar for the three 
samples. Its profile differs in the two individual samples 
(FWHM(NLSls) = 2050 km s" 1 , FWHM(z > 4 QSOs) = 
3100 km s _1 ); as might be expected, the NLS1 core contri- 
bution is narrower than that in the high z QSOs. The sec- 
ond component, PC2, accounts for spectrum-to-spectrum 

n reported in several previous studies; see Osmer & Shields (1999) for 
comparison. 
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variations present in the wings of the lines. This is again 
exhibited by both samples but in a much larger proportion 
by the NLSls. Also, the line-core and the line-wing modu- 
lations show opposite trends, suggesting that objects with 
prominent wings possess weak cores. Higher-order com- 
ponents continue to display, in different proportions, the 
core and the wing variations, but their intricacies become 
difficult to characterize and compare. 

Associating physical explanations to the properties ex- 
hibited by PCs is difficult, but additional information can 
be gained from the statistics of the PCA results. Table 
4 lists the proportions (fractions of total variance) con- 
tributed by the first five PCs for each individual sample 
and for the combined one. The numbers indicate that the 
NLS1 sample can be better represented than the z > 4 
quasar sample by a low number of principal eigenvectors. 
The difference in the PCs' individual and cumulative pro- 
portions indicate that NLS1 sources comprise a more spec- 
troscopically compact sample than the high z quasars. 
Another basis for evaluating these results is a statistical 
comparison of the distributions of the weights of the first 
two principal components for the NLS1 and z > 4 QSO 
objects. In the PCA run for the combined sample, i.e., 
for a common set of eigenvectors representing the NLSls 
and the high z QSOs, the Kolmogorov-Smirnov (KS) test 
shows that the possibility of the two data sets being drawn 
from the same parent population is not excluded, but the 
evidence for this is weak: the KS probabilities are 0.047 
and 0.188 for weights of PCI and PC2 respectively. 

As a consistency test, we also performed this analysis 
for a much larger spectral range, covering Lya to He II 
bandwidth, but with a lower number of NLSls (only 16 
objects), and therefore lower statistics; the main results 
remained unchanged. NLSls and z > 4 QSOs are proba- 
bly not close spectroscopic or physical analogues. 

5. CONCLUSIONS 

In this paper we present an analysis of all publicly avail- 
able spectra for NLS1 galaxies in the HST archive. The 
resulting sample of 22 NLSls with spectra spanning the 
UV-blue wavelength range is the largest that has been 
used, to date, in emission line studies of these objects. 
We employed these data to construct composite spectra 
(average and median), and thus, to characterize in detail 
the typical spectral properties of the NLS1 class. 

The resulting composite spectra are used to estimate the 
strengths of a large number of emission lines, and to quan- 
tify the continuum shape over a broad bandpass. Power- 
law fits to the continuum indicate a discrepancy from the 
results obtained from more general AGN composite stud- 
ies: NLSls have steeper UV-bluc spectra. Possible ex- 
planations include intrinsic reddening and a trend of the 
low redshift sources to have intrinsically softer continua. 
A relation between the continuum shape and the redshift 
is not readily evident in this sample. A significant cor- 
relation is however observed between the spectral slope 
and luminosity, indicating that the redness of the NLSls 
is related to the low luminosity of these objects. More- 



over, the apparent connection between the UV resonance 
absorption lines and luminosity suggests that the steep 
slopes measured in these objects are due at least partially 
to reddening. The ionization state of the absorbing ma- 
terial and its relationship to the accretion source are not 
well determined, however. 

The NLS1 composites additionally allow us to quantify 
emission-line velocity offsets. The correlation between the 
velocity shifts and the degree of ionization that is found 
in normal broad- line AGNs is also present in the NLS1 
sources, in both permitted and forbidden lines. This re- 
sult may be of interest for comparison with NLS1 model 
predictions. 

The NLS1 data permit further investigation into the 
proposed analogy between these sources and the z > 4 
quasars. Previous work suggested that both may be de- 
scribed by a high accretion rate and super-solar metallic- 
ities. The comparative study that we conduct based on 
the NLS1 and high z quasar average UV emission prop- 
erties reveals a significant contrast between their spectral 
characteristics. The composite spectra exhibit primarily 
the anticipated differences associated with the Baldwin 
Effect, i.e., relatively stronger lines in the less luminous 
sources. Nonetheless, the comparison makes it evident 
that the metal enrichment of the surrounding gas in the 
high z QSOs is higher than that in the NLS1 objects. This 
result does not support the hypothesis that these two types 
of objects are similar in their detailed physical character- 
istics or evolutionary phase. 

Additional confirmation of the differences exhibited by 
the two types of objects comes from a principal component 
decomposition applied directly to the spectra. The spec- 
tral PC analysis indicates that the z > 4 sources are sta- 
tistically likely to share only part of their emission proper- 
ties with the NLS1 galaxies; in particular, the low-velocity 
component, although enhanced in the high z QSOs, is less 
prominent than in the NLS1 objects. The statistics show 
that the z > 4 QSO phenomenon is controlled by a broader, 
more heterogeneous family of properties than those gov- 
erning the NLSls; the NLS1 spectroscopic features can be 
reconstructed from a smaller number of parameters (eigen- 
vectors) than that necessary in the case of z > 4 sources. 
In conclusion, NLSls and high z quasars are spectroscop- 
ically disparate, and it is therefore doubtful that a close 
physical connection exists between these source types. 
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Fig. 1. — Histogram of continuum spectral indices, a, where F v oc v a , for all but one (MrkllO) NLS1 in the sample. The average and 
median values arc indicated. 
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Fig. 2. — Top panel: NLS1 composite spectrum plotted as log(F^) vs. rest-frame wavelength, with the principal emission features identified. 
The flux has been normalized to unit mean flux over the wavelength range 1430 A- 1470 A. Middle panel: The median composite spectrum, 
plotted on a linear scale, and zoomed near the continuum level for a better visualization of the weak features in the optical range. A more 
detailed UV line identification is presented in Figure 4. The power-law continuum fits are overplotted as dashed and dotted lines. Bottom 
panel: Number of NLSls contributing to the composite as a function of rest-frame wavelength. 
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Fig. 3. — Spectral indices plotted vs. redshift (upper panel) and 1450A luminosity (lower panel), with L v expressed in ergs s — 1 Hz -1 . The 
Spearman rank coefficient and the probability of the correlation happening by chance are indicated. The error bars in both directions are 
smaller than the symbol size, and therefore not indicated. 
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Fig. 4. — Top panel: NLS1 average composite spectrum plotted for the UV range only, where at least 18 objects are contributing. The same 
normalization as in Figure 2 is used. The prominent emission features are marked. Middle panel: Flux standard deviation in F^ relative to 
the average composite spectrum, as a function of rest-frame wavelength for the individual spectra comprising the composite. Bottom panel: 
The standard deviation of the mean, expressed in percent. 
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Fig. 5. — The optical range of the HST NLS1 median composite, constructed using only the 3 objects that span the whole spectral range. 
The same flux normalization as in Figure 2 is used. The Sulentic et al. (2002) NLS1 median composite, based on 24 ground-based object 
spectra, is shown for comparison. The strong similarity between the two medians suggests that the ffiJTNLSl composites are representative 
of these objects. 
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Fig. 6. — Emission-line velocity offsets, relative to the rest frame (defined by C IVA1549), as a function of ionization potential for selected 
emission lines (see text). Error bars show the 1 <r uncertainty in the velocity measurement. Permitted and semi-forbidden lines are shown in 
the top panel, and forbidden lines are shown in the bottom panel. Only one measurement is shown for each line, and the point type indicates 
the adopted profile. Overlapping points are slightly offset horizontally from each other for clarity. The points are labeled by ion. 
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Fig. 7. — Direct comparison between NLS1 and z > 4 QSO average composite spectra. The same normalization as in Figure 2 is used. 



17 



3.5 



3.0 - 



-> 2.5 h 
o 

$ 2.0 - 

Ld 
cn 

o 1.5 h 



1.0 h 
0.5 



o o 
o 



oNLSIs 

• z > 4 QSOs 



H = 50, q = 0.5, A = 



o 
o 



"n-O. 



-0 



• w • 
• • • # 



26 



28 30 
log L u (1450 A) 



32 



Fig. 8. — Rest-frame EW of the C IV emission line, in A, as a function of 1450A luminosity, in ergs s _1 Hz -1 , for individual NLS1 and 
2 ^ 4 QSO spectra. The Baldwin relation found by Osmer, Porter & Green (1994) for a sample of 186 luminous quasars [logL„(1450) > 29, 
< z < 3.8] is shown as a dotted line. 
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Fig. 9. — Mean and standard deviation spectrum and the first 5 principal components given by PCA, performed for the C IV region 
only, applied to the NLS1 sample (dotted line), the z > 4 QSOs (dashed line), and the combined sample (continuous line). The principal 
components show similar modulations in the spectral variation in both categories of objects; the amplitude differs. 
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Table 1 
NLS1 Sample 



Object 


Instr. 


Gratings 


Coverage 3, 


Ref b 


Ark 564 


FOS 


G130H,G190H,G270H,G400H,G570H 


1087-6817 


6,9,10,11 


1H0707-495 


STIS 


G140L,G230L 


1122-3155 


6 


IRAS 13224-3809 


STIS 


G140L.G230L 


1117-3148 


1,6 


IRAS13349+2438 


FOS 


G190H,G270H 


1572-3294 


4 


KUG 1031+398 


FOS 


G130H,G190H,G270H 


1087-3301 


8,11 


MrkllO 


STIS 


G140M 


1194-1250 


2,11 


Mrk335 


FOS 


G130H,G190H,G270H 


1087-3301 


2,11 




GHRS 


G160M 


1221-1257 




Mrk478 


FOS 


G130H,G190H,G270H 


1087-3277 


2,5,11 




STIS 


G140M 


1194-1300 




Mrk486 


FOS 


G190H,G270H 


1567-3293 


2,11 


Mrk493 


FOS 


G130H,G190H,G270H,G400H,G570H 


1087-6817 


7,11 


Mrk766 


STIS 


G140L,G230L 


1087-3151 


3,5,7,11 


NGC4051 


STIS 


E140M 


1140-1729 


6,11 


PG 1211+143 


FOS 


G130H,G190H,G270H 


1087-3275 


2,11 




GHRS 


G140L 


1190-1477 






STIS 


G140M 


1194-1300 




PG 1351+640 


FOS 


G130H,G190H,G270H 


1087-3301 


2 




STIS 


G140M 


1194-1300 






QTTQ 

alio 




lODo-olOl 




PG 1404+226 


FOS 


G130H,G190H,G270H,G400H 


1087-4780 


2,11 


PG 1411+442 


FOS 


G130H,G190H,G270H 


1087-3276 


2 




STIS 


G230L 


1572-3157 




PG 1444+407 


FOS 


G130H,G190H 


1087-2330 


3 


RX J0134-42 


FOS 


G130H,G190H,G270H,G400H,G570H 


1087-6817 


5 


Ton S180 


STIS 


G140M 


1194-1299 


5,11,12 




STIS 


G140L,G230L 


1120-3160 




WPVS007 


FOS 


G130H,G190H,G270H,G400H,G570H 


1087-6817 


5,13 


I Zwl 


FOS 


G130H,G190H,G270H 


1087-3276 


2,11 




FOS 


G190H,G270H 


1568-3295 






FOS 


G270H 


2222-3277 






GHRS 


G160M 


1221-1257 




II Zwl36 


STIS 


G140M 


1194-1300 


2 




GHRS 


G140L 


1153-1739 





a given in A 

b References to optical emission-line measurements that led to their NLS1 classification 



References. — (1) Boiler et al. (1993); (2) Boroson & Green (1992); (3) Goodrich (1989); 
(4) Grupe et al. (1998); (5) Grupe et al. (1999); (6) Leighly (1999); (7) Osterbrock & Pogge 
(1985); (8) Puchnarewicz et al. (1995); (9) Stirpe et al. (1990); (10) van Groningen (1993); 
(11) Veron-Cetty, Veron & Goncalves (2001); (12) Winkler (1992); (13) Winkler et al. (1992) 
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Table 2 



NLS1 Properties 


Object 


z a 


log J L„(1450) b 


a(f„ oc v a ) 


EW(C IV) C 


Trosat 


Ref 3 


ATK 004 


U.uzoy 


9£ 9^ 


1 Q/\_i_n, ni 
- 1 .04xU.Ul 


ol .0 


o.4ztU. 1 


i 
i 


InlJ i U I -4y0 


U.UoOO 


90 ni 


-U.40±U.U1 


9ft £ 


Z.O±U.O 


Q 
O 


TR Am "?994 ^XflQ 


u.uoou 


9Q 


44+0 09 





4 p;-i-n 1 

4.0ICU. 1 


1, O 


TO A«1 '?'?4q-L94'?8 


U.1UOO 


9R ^ 






o e-i-n i 

z . oznu . i 


Q 
O 




fl 04 "}9 


98 10 




Q9 4 

fyZ .4: 


4 1 


1 

J. 


1V1I K11U 


U.UulU 








Z.4znU. 1 


/I 

4 


Mrk335 


0.0260 


29.23 


-0.64+0.01 


77.7 


2.9+0.1 


3 


Mrk478 


0.0755 


29.86 


-0.68+0.01 


39.7 


3.1+0.1 


1, 3 


Mrk486 d 


0.0387 


27.84 


-2.84+0.09 


50.0 






Mrk493 


0.0311 


28.35 


-0.73+0.01 


107.5 


2.7+0.2 


1 


Mrk766 d 


0.0127 


26.13 


-4.41+0.04 


98.1 


2.7+0.1 


1, 3 


NGC4051 d 


0.00245 


26.38 


-1.86+0.04 


75.5 


2.8+0.0 


3 


PG 1211+143 


0.0804 


29.75 


-1.14+0.01 


66.5 


3.1+0.2 


3,4 


PG 1351+640 d 


0.0875 


29.65 


-1.35+0.01 


115.0 


2.5+0.6 


4 


PG 1404+226 d 


0.0935 


29.22 


-0.91+0.01 


44.7 


4.1+0.2 


3,4 


PG 1411+442 d 


0.0875 


29.68 


-0.98+0.03 


148.0 


3.0+0.5 


4 


PG 1444+407 


0.2640 


30.55 


-0.81+0.04 


42.1 






RX J0134-42 d 


0.2360 


30.31 


-0.12+0.01 


27.5 


7.7+2.6 


2 


Ton S180 


0.0613 


29.63 


-0.76+0.01 


33.0 


3.0+0.1 


3 


WPVS007 d 


0.0280 


27.78 


-1.41+0.01 


112.2 


9.0+2.0 


2 


I Zwl 


0.0585 


29.40 


-1.75+0.01 


50.0 


3.1+0.1 


1, 3 


II Zwl36 


0.0619 


29.59 


-0.73+0.04 


61.2 


3.2+.. . 


4 



a calculated from the observed spectra using the C IV emission-line; the only exception is 
MrkllO for which the line is not available, and the redshift is taken from the literature 

b to ease comparison with earlier published work, an H = 50 km s _1 Mpc _1 , q = 0.5, 
A = cosmology was adopted. The values reflect flux measurements which are corrected for 
Galactic extinction; typical uncertainties are ^0.05 dex. 

c rest-frame EWs (in A) calculated by fitting a single Lorentzian profile to the line 

d Emission lines are contaminated by absorption; indicated measurements (EWs), obtained 
from the polynomial interpolated profiles, should be regarded as lower limits 

e References for the Trosat values 



References. - (1) Boiler, Brandt & Fink (1996); (2) Grupe et al. (1998); (3) Leighly 
(1999); (4) Wang et al. (1996) 
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Table 3 



Emission-line measurements 





A(ab' 

(A) 


Z 


X , d 

^mean_rest 


Aw 


Rel. Flux f 


EW g 


FWHM h 


Line a 




(A) 




(km s- 1 ) 


[100 xF/F(Lya)l 


(A) 


(km s- 1 ) 


C III* 


1175. 


.5 


1175 


.5±0 


.0 


0±130 


0.9±0.2 


1.1 


1190± 690 


Si II 


1197. 


,4 


1199 


.3±0 


7 


480±170 


0.5±0.2 


0.6 


790± 500 


Si III 


1206. 


,5 


1206 


.7±1 


Q 


50±470 


0.4±0.6 


0.5 


780±1350 


Lya 


1215. 


.7 










100.0±1.6 


131.8 




..narrow 


1215 


.7 


1216 


.3±0. 


i 

. j. 


150±25 


69.5±0.9 


91.6 


1480± 10 


..broad 


1215 


.7 


1212 


,6±1 


A 


-760±340 


30.6±1.0 


40.2 


7000± 280 


N V 


1240. 


.8 


1240 


.0±0 


1 

. 1 


-190±20 


13.0±1.5 


17.6 


2110± 130 


Si II* 


1248. 


.4 


1245 


.0±0 


A 

. -i 


-810±90 


3.7±1.0 


5.0 


2200± 400 


Si II 


1264. 


.8 


1262 


.9±0. 


Q 

. iJ 


-450±210 


3.5±0.4 


4.8 


3510± 350 


01 


1303 


.5 


1305 


,4±0. 


Q 
.0 


440±70 


3.7±0.6 


5.5 


2060± 130 


Si II 


1309. 


.3 


1309 


,1±1. 


r. 

• O 


-40±340 


1.0±0.7 


1.5 


2030± 280 


CII 


1335 


.3 


1336 


,5±0. 


Q 
.O 


270±70 


2.6±0.3 


4.0 


1930± 260 


N II b 


1344. 


,6 


1349 


,0±0. 


O 

,y 


980±200 


1.1±0.3 


1.7 


2020± 430 


Si IV+O IV] 


1400 


.0 










10.5±3.2 


17.2 




..Si IV 


1393 


.7 


1393 


,7±0 


7 


0±150 


3.1±2.2 


5.1 


1520± 600 


..0 IV] 


1401 


.4 


1404 


,6±1. 


,0 


680±320 


5.4±2.3 


8.9 


2030± 440 


..Si IV 


1402. 


.7 


1399 


.9±2. 


7 


-590±570 


2.0±0.7 


3.3 


1520± 690 


N IV] 


1486 


.5 


1486. 


.8±0. 


Q 

.y 


60±180 


0.7±0.3 


1.2 


2020±1910 


Si II 


1530 


.1 


1536 


,9±0. 


.Z 


1330±40 


4.6±0.8 


7.9 


2030± 280 


C IV 


1549. 


,5 










44.8±1.7 


77.1 




..narrow 


1549 


.5 


1549 


,5±0 


1 

.1 


0±20 


33.3±1.3 


57.3 


1900± 60 


..broad 


1549 


.5 


1541 


,6±8. 


,z 


-1520±1580 


11.5±1.1 


19.8 


12000±2500 


He II 


1640. 


.4 










23.3±2.6 


41.5 




..narrow 


1640 


.4 


1639 


,7±0 


.z 


-120±40 


6.4±0.8 


11.4 


1780± 240 


..broad 


1640 


.4 


1630 


.4±3 


n 
u 


-1820±540 


16.9±2.5 


30.1 


11990±4090 


III] 


1663 


.5 


1665 


,4±0 


.z 


340±40 


4.1±1.1 


6.2 


2500± 480 


Al II 


1670 


.8 


1680 


.0±1 


c; 
o 


1650±260 


1.3±1.0 


1.9 


2250±1330 


Nil 


1725. 


.2 


1725 


,6±0 


,y 


70±150 


2.2±0.9 


3.3 


4710±2100 


N III] 


1750 


.5 


1750. 


.5±0 


.6 


0±50 


2.0±0.2 


3.0 


1480± 170 


Fe II UV191 


1785. 


.4 


1785 


.3±0 


7 


-20±110 


1.2±0.2 


1.8 


1430± 250 


Si II 


1814. 


.1 


1816 


,9±0. 


.0 


460±80 


0.6±0.2 


0.9 


1250± 730 


Al III 


1858. 


.8 


1853 


,5±1. 


1 

. J. 


-850±170 


1.3±0.5 


1.9 


2240± 540 


Fe III UV52 


1867. 


.9 


1864. 


,1±0 


Q 
. O 


-610±120 


1.7±0.5 


2.6 


2760± 730 


Si III] 


1892 


.0 


1892. 


,7±0. 


1 
. 1 


110±20 


3.6±0.5 


5.5 


1250± 130 


C III 


1907. 


.9 










21.1±5.7 


32.2 




..narrow 


1907 


.9 


1908 


,2±0 


1 

. 1 


50±20 


10.0±0.6 


15.3 


1380± 10 


..broad 


1907 


.9 


1907 


.7±4. 


. u 


-30±720 


11.0±5.7 


16.9 


12230±1110 


N II] 


2141. 


.4 


2144 


.0±0 


Q 


360±130 


0.3±0.2 


0.4 


1040±1140 


CII] 


2327 


.5 


2326 


.6±0 


Q 
-O 


-120±100 


3.7±1.3 


5.7 


2540±430 


Fe III UV 47 


2419. 


,3 


2422 


.8±0 


9 


430±30 


3.2±0.5 


5.1 


1560±200 


on 


2438. 


.8 


2439 


.3±0 


.0 


60±70 


0.7±0.2 


1.2 


1060±250 


[0 II] 


2470. 


,9 


2467 


,6±0 


.8 


-400±100 


0.9±0.3 


1.5 


1570±470 


Cl] b 


2478. 


.6 


2484. 


.0±0 


.1 


650±10 


2.3±0.6 


3.7 


2580±510 


Al II] 


2669. 


.9 


2674 


,7±0 


.4 


540±50 


0.8±0.6 


1.3 


1050±430 


Mg II 


2797. 


.9 










31.4±6.6 


58.0 




..narrow 


2797 


.9 


2801 


,6±0. 


.3 


390±30 


11.2±2.2 


20.7 


1680±60 


..broad 


2797 


.9 


2798 


,3±0 


.4 


40±40 


20.2±6.5 


37.4 


4620±1090 


III 6 


2960. 


,6 


2957 


,4±1. 


.7 


-320±170 


0.3±0.2 


0.9 


850±600 


Fe IP 


2964. 


.3 


2966 


,7±3 


.6 


240±360 


0.1±0.1 


0.2 


830±900 


IV° 


2982. 


.5 


2981 


.9±2 


.5 


-60±250 


0.3±0.3 


0.9 


820±610 


Ne III] 


2986 


.9 


2990 


.0±0 


.1 


310±10 


0.2±0.1 


0.5 


830±680 


F V° 


3109. 


,0 


3109 


,5±0. 


.4 


50±40 


0.6±0.1 


2.0 


550±10 
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Table 3 — Continued 





Mab 


z 


\ d 

^mean^est 


A 

l\v 


Kel. l^mx 




Till TTT1\ Til 

b WH1V1 


Line a 


(A) 




(A) 




(km s ) 


[100 xF/F(Lya)] 


(A) 


(km s ) 


IIP 


3122 


.5 


3122 


,9±1 


Q 
. o 


40±120 


0.8±0.4 


2.7 


910±370 

\J _i_ \J — [ — U 1 u 


IIP 


3133 


.7 


3134 


,2±1 


Q 
. o 


50±120 


1.8±0.4 


6.0 


940±210 


C IP 


3166 


.7 


3165 


.3±0 


.6 


-130±60 


0.8±0.3 


2.7 


760±390 


He I 


3188. 


.6 


3190 


,6±1 


,7 


190±160 


0.7±0.1 


1.6 


990±230 


IV b 


3216 


.8 


3215 


.7±0 


.5 


-100±50 


0.5±0.1 


1.1 


830±120 


Ne II b 


3230 


.6 


3230 


,9±0 


s 

.O 


30±70 


0.9±0.2 


2.1 


940±260 


Fe I Opt91 


3261 


.2 


3260 


,8±0. 


7 
. f 


-40±60 


1.0±0.1 


2.5 


840±130 


Fe II Optl 


3281 


.2 


3281 


,6±0 


.5 


40±40 


1.2±0.1 


3.0 


840±70 


[Fe III] b 


3308 


.5 


3307 


,8±2 


A 

.4: 


-60±220 


0.7±0.1 


1.9 
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Table 3 — Continued 







\ d 


Av 


Rcl. Flux f 


EW g 


FWHM h 


Linc a 


(A) 


(A) 


(km s -1 ) 


[100 xF/F(Lya)] 


(A) 


(km s _1 ) 


[FeX] 


6376.3 


6373.5±0.6 


-130±30 


1.2±0.1 


6.5 


810±90 


[N II] 


6548.8 


6549.5±0.9 


30±40 


1.8±0.6 


10.1 


760±70 


Ha 


6564.6 






61.0±1.4 


344.2 




..narrow 


6564.6 


6565.9±0.1 


60±5 


57.8±1.0 


325.8 


760±20 


..broad 


6564.6 


6566.7±0.4 


90±20 


3.3±1.2 


18.4 


5500±1140 


[N II] 


6585.3 


6585.9±0.9 


30±40 


6.1±0.6 


34.3 


760±70 



a Emission lines are measured from the average composite spectrum unless otherwise 
noted 

b Uncertain identifications 

c vacuum wavelengths, both below and above 2200A 
d given by the centroid of the Lorentzian fit 

'The line fluxes are normalized to the Lya flux to ease comparison with published 
results for other AGN/quasar samples 

c Lines measured from the median composite 

s The errors in EWs are not quoted as the relative errors are the same as for the fluxes 

h Intrinsic width of the emission line; an instrumental broadening of 230 km s _1 , the 
approximate FOS spectral resolution, high-resolution gratings, was assumed 



Table 4 

Proportions of PCs in the C IV emission region 



Component 


NLS1 


z > 4 QSO 


Combined sample 




(21 obj.) 


(44 obj.) 


(65 obj.) 


PCI 


0.358 


0.272 


0.296 


PC2 


0.139 


0.077 


0.091 


PC3 


0.097 


0.071 


0.061 


PC4 


0.064 


0.053 


0.051 


PC5 


0.046 


0.044 


0.041 



